Magnetic anisotropy effects on quantum impurities in superconducting host 
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We study the magnetic anisotropy effects on the localized sub-gap excitations induced by quantum impurities 
coupled to a superconducting host. We establish the ground- state phase diagrams for single-channel and two- 
channel high- spin Kondo impurities; they unveil surprising complexity that results from the (multi-stage) Kondo 
screening in competition with the superconducting correlations and the magnetic anisotropy splitting of the spin 
multiplets. We discuss the possibility of detecting the Zeeman splitting of the sub-gap states, which would 
provide an interesting spectroscopic tool for studying the magnetism on the single-atom level. We also study the 
problem of two impurities coupled by the Heisenberg exchange interaction, and we follow the evolution of the 
sub-gap states for both antiferromagnetic and ferromagnetic coupling. For sufficiently strong antiferromagnetic 
coupling, the impurities bind into a singlet state that is non-magnetic, thus the sub-gap states move to the edge of 
the gap and can no longer be discerned. For ferromagnetic coupling, some excited states remain present inside 
the gap. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.20.Hb, 73.20.-r 



I. INTRODUCTION 

Tunneling spectroscopy is the prevalent experimental ap- 
proach for studying superconductivity. It provides informa- 
tion on such fundamental properties as the energy gap, pair- 
ing symmetry, and pairing interactions 1-6 . Using a scanning 
tunneling microscope (STM) it is possible to examine the im- 
purity effects on the single-atom level 4 ^ - — ; see Fig.Hk). Such 
measurements provide crucial data on the nature of the super- 
conducting state in complex materials 3,5 . With improvements 
in the instrumentation, experiments are being performed at in- 
creasingly low temperatures and ever further details in the lo- 
cal density of states (LDOS) can be resolved: recent STM 
work performed in the 300 mK range on magnetic adatoms 
adsorbed on superconductors has clearly revealed the exis- 
tence of multiple sub-gap excitation peaks^r— ; see Fig.[TJ)). 
It was proposed that these may be interpreted as the magnetic- 
impurity-induced bound states associated with the different 
angular-momentum scattering channel s 19 i 22 i 23 , but we show 
in this work that an alternative interpretation in terms of the 
magnetic-anisotropy effects is also possible. 
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Figure 1 : (Color online) (a) Magnetic impurity atom adsorbed on the 
superconductor surface probed by the STM. (b) Idealization of the 
characteristic differential conductance spectrum as recorded by the 
STM above a magnetic impurity: in the gap, the LDOS is essentially 
zero except for the discrete peaks which correspond to the transitions 
from the ground state to the sub-gap excited states. 



Since the magnitude of the magnetic anisotropy of mag- 
netic adsorbates is comparable or may even exceed the su- 
perconducting gap 24,25 , its effects are drastic, yet they have 
received hardly any attention in this context. In addition, it 
has been very recently demonstrated that the strength of the 
anisotropy (parameter D) in magnetic molecules can be con- 
tinuously tuned in mechanical break junctions 26 ; using super- 
conducting leads, one could then directly study the effect of 
the magnetic anisotropy on the transport and the spectral prop- 
erties of magnetic molecules coupled to superconductors. Fi- 
nally, spin-orbit coupling is always present for any impurity 
atom embedded in the superconducting bulk, therefore the 
class of problems in which an anisotropic spin degree of free- 
dom is coupled to a superconducting host is indeed wide rang- 
ing. 

The theory of magnetic impurities in a superconductor was 
first worked out within simplified classical-spin model a 27 i 28 , 
while later works took into account the quantum nature of 
impurities and, among other improvements, properly de- 
scribed the competition between the screening of the impu- 
rity moment by the Kondo effect and the superconducting 
correlations 29 . The crucial difference between a classical and 
quantum spin is that a classical spin is static (essentially equiv- 
alent to a local point-like magnetic field), it has no internal 
dynamics, and it cannot flip the spin of the conduction-band 
electrons. For this reason, the impurity problem in a classical- 
spin approximation is a non-interacting (mean-field) theory. 
The quantum impurity, however, needs to be treated using 
the tools of the many-particle theory, which can take into ac- 
count the non-perturbative effects induced by strong interac- 
tions. The behavior of classical and quantum impurities is 
very different in many respects. For classical impurities, there 
is no difference if the coupling to the conduction-band elec- 
trons is ferromagnetic or antiferromagnetic, while real quan- 
tum impurities have very different properties in these two 
cases. Furthermore, in the absence of magnetic anisotropy, 
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a quantum impurity model has full SU(2) spin symmetry and 
any breaking of this rotational invariance would indicate a de- 
ficiency of the method used; a classical spin, however, explic- 
itly breaks the spin symmetry even at the level of the Hamilto- 
nian itself. Unless there is a physical mechanism that can lead 
to a real symmetry breaking [such as spontaneous symmetry 
breaking due to magnetic ordering induced by inter-impurity 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction], only a 
quantum impurity model can thus provide the qualitatively 
correct result, while a classical model will be affected by 
spurious symmetry breaking. (In magnetically-ordered sys- 
tems, classical spin models may be fully adequate for many 
purposes.) The spin SU(2) symmetry breaking also implies 
the breaking of the time-reversal symmetry. This has impor- 
tant consequences for the degeneracy of the sub-gap states 
(cf. Kramers' degeneracy theorem) and it is one of the ma- 
jor differences between quantum and classical (static) impu- 
rities. This dichotomy recently gained renewed attention in 
the context of topological insulators^"—, since only impuri- 
ties which break the time-reversal invariance can open the gap 
in the topologically-protected edge/surface states^ - —. 



II. MODEL AND METHOD 

We describe the impurity system by the Hamiltonian 
H = i^band + i^imp, where i^band describes N channels of 
conduction-band electrons using the mean-field BCS Hamil- 
tonian with the gap 
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while the impurity is described by a Kondo-like Hamiltonian 
with a magnetic anisotropy terrr 4 ^^- - — • 
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Here Ji are the exchange coupling constants, Si mp is 
the impurity spin operator satisfying the su(2) Lie algebra 
[<Sim P ,a, Si mPi p] = ie a p<ySi mpn , and s$ is the channel-i spin 
density at the position of the impurity: 



Accurate calculations for quantum impurities in supercon- 
ductors became possible by generalizing the numerical renor- 
malization group (NRG) method to problems with a supercon- 
ducting electron batr42. This method was applied to simpli- 
fied models such as the spin- 1/2 Kondo modeler— and the 
non-degenerate Anderson impurity model 46 . Real magnetic 
impurity atoms and molecules generally require, however, a 
more sophisticated description in terms of the multi-channel 
degenerate Anderson model or the high- spin Kondo model 
with magnetic anisotropy terms 4 ^. The goal of the present 
work is thus to apply the NRG to study the sub-gap excitations 
for such more realistic models, featuring anisotropic high- 
spin Kondo impurities with more than one channel. There 
are some related works in the literature. Lee et al. have 
studied the isotropic high- spin impurities in a side-coupled 
configuration 48 ' and the two-level impurity model 49 . The two- 
channel models have been studied in the context of unconven- 
tional superconductors 50 and iron pnictide superconductors^. 
Moca et al. have demonstrated the presence of multiple 
sub-gap states in a multi-orbital model for Mn impurity in 
MgB2 52 . Multiple sub-gap states also appear in the case of 
XXZ anisotropic Kondo exchange coupling in the ^imp = 1/2 
Kondo model 4 ^. Multi-channel high- spin Kondo models with 
magnetic anisotropy terms, however, have not yet been stud- 
ied. 

The paper is structured as follows. In Sec. HH we define 
the model and comment on its relevance for actual adsorbate 
systems. In Sec. |III1 we study the ground state diagrams for 
various one- and two-channel isotropic high- spin Kondo sys- 
tems, while the magnetic anisotropy effects are presented in 
Sec. [IV] The spectral peaks and their splitting due to the 
magnetic anisotropy in the multi-channel case are discussed in 
Sec.[V] Section|Vl]is devoted to the role of the external mag- 
netic field, which splits all many-particle states with S ^ 0. 
Finally, impurity dimers are studied in Sec. I VIII 
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where M is the number of the states in the conduction band. 
Furthermore, D is the longitudinal anisotropy, g is the im- 
purity ^-factor, fiB is the Bohr magneton, and B is the ex- 
ternal magnetic field (non-zero B is discussed in Sec. [VTb . 
The multiple channels correspond to the different symmetry 
channels of the Bloch states, which hybridize with the impu- 
rity d-levels; our high- spin Kondo model may be thought to 
arise from some multi-orbital Anderson model after perform- 
ing the Schrieffer- Wolff transformation^"— . Strictly speak- 
ing, the low-energy effective model for a multi-orbital Ander- 
son model is not necessarily a high-spin Kondo model. Such 
an exception occurs, for example, if the system is not in the 
local-moment regime, but rather its valency is strongly fluctu- 
ating, or if there is also some orbital moment on the impurity 
atom. For surface- adsorbed impurities, the symmetry in real 
space is broken, thus it is reasonable to expect strong quench- 
ing of the orbital moment. The valence- fluctuation regime 
cannot be excluded a priori, but the systems of such complex- 
ity are beyond the capabilities of the NRG method. In this 
work we thus focus on the problems where the orbital moment 
is quenched and the electrons in the d-orbitals are locked into 
a high- spin state by the strong Hund's coupling. Such cases 
are adequately described by the proposed model. For D ^ 
or B ^ 0, the Hamiltonian only has an axial U(l) spin sym- 
metry, thus the sole conserved quantum number is S z , the z 
component of the total spin. 

The NRG method consists of discretizing the continua of 
the conduction-band states, rewriting the Hamiltonian in the 
form of one-dimensional tight-binding chains with an expo- 
nentially decreasing hopping constants, and diagonalizing the 
resulting Hamiltonian iteratively by adding one further chain 
site per channel in each step£&~— . The spectrum of many- 
particle states is truncated to the low-energy part after each 
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step. Due to the low symmetry of the problem (there is no 
particle conservation in the superconducting case, and there 
is only partial spin symmetry in the presence of magnetic 
anisotropy or magnetic field), these calculations are numeri- 
cally very demanding. The size of the matrices that need to 
be diagonalized at the given truncation cutoff strongly depend 
on the value of the discretization parameter A. In many situa- 
tions, one can use a large value of A and reduce the discretiza- 
tion artifacts by the so-called z-averaging trickSr— ; this ap- 
proach produces excellent results for featureless (flat) conduc- 
tion bands. One needs to be very careful, however, in the 
vicinity of phase transitions, since the ground state obtained 
in a calculation can be z-dependent (in other words, it is pos- 
sible that for exactly the same model parameters, one obtains 
a different ground state for different interleaved discretization 
meshes). The averaging may then be ill-defined for an interval 
of parameters where such z-dependence of the ground state 
occurs. The width of this interval grows as A is increased. 
Nevertheless, experience shows that in spite of this difficulty, 
one can determine the transition point very accurately even by 
performing the calculations with a very large value of A. Test 
calculations on simple problems, for example, show that by 
determining the value of the system parameter where the tran- 
sition occurs for a fixed value of z, and averaging such results 
over z, one obtains a value that changes little with A. Often 
it is sufficient to use only two values of z (such as 1 and 0.5) 
to obtain good results. This approach has been used, for ex- 
ample, to establish the accurate values tabulated in Table U in 
the following section. In other parts of this work, where high 
accuracy was not essential, we performed no such averaging, 
thus the results are only qualitatively correct. 



III. TRANSITIONS IN THE ISOTROPIC MODEL 

The properties of the multi-channel Kondo model in the 
normal state depend on the relation between the number of 
channels and the impurity spin; roughly speaking, each chan- 
nel can screen one half unit of the impurity spinal - —. Thus, 
for TV < 25'i m p, the impurity spin can be only partially 
screened, while for TV = 2S{ mp there is an exact spin compen- 
sation, yielding a singlet ground state (GS); for TV > 2Si mp , 
exotic non-Fermi-liquid (NFL) states may arise^. In the 
superconducting state, the behavior of the Kondo model is 
only well explored for the simplest case of TV = 1 and 
Simp = l/2fi: as the gap A is increased from zero, there 
is a transition at A c between the regime where the impurity 
is Kondo- screened (singlet ground state) to the regime where 
the impurity is free (doublet ground state). The value of A c 
is of the order Tk, the Kondo temperature, which is the char- 
acteristic energy scale of the Kondo effect. The transition oc- 
curs because in the superconducting state there is an insuffi- 
cient number of low-energy electron states to participate in the 
Kondo screening. 

We now consider the general multi-channel high-spin case, 
first in the absence of the anisotropy and for equal coupling 
constants Ji for all channels. By analogy with the known re- 
sults, we expect that for 26 f i mp > TV there is a transition from 



the Kondo screened S = 5i mp — TV/2 state to the "free-spin" 
S = S'imp state as the gap is increased. Our NRG calculations 
fully support this picture. For 2S[ mp < TV the NFL effects 
make a-priori predictions difficult; in numerical simulations 
for TV = 2 and 5i mp = 1/2, we observe a transition from a 
degenerate pair of singlet states to a doublet ground state as A 
increases. The A c for TV = 1 and TV = 2 are tabulated in Ta- 
ble U for a range of Si mp . The ratio A c /Tk strongly depends 
on Si mp in spite of the fact that the superexchange couplings 
are constant and thus the Kondo scale is formally the same 
in all cases. The variation with TV is weaker, with a notable 
exception of the NFL case with TV = 2, = 1/2, where 
A c /Tk is much reduced. 
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Table I: Dependence of the ratio A c /Tk between the critical gap 
and the Kondo temperature on the number of channels TV and the 
impurity spin ftmp. In the two-channel case, the coupling to both 
channels is taken to be equal, Ji = J. In all calculations J — 0.2, 
and the Kondo temperature is Tk = 1.16 x 10 -5 . We use Wilson's 
definition of the Kondo temperature^ throughout this work. 

It must be emphasized that even in the limit of small ex- 
change couplings that is, when the Kondo effect plays no 
role and the impurity remains unscreened, the system still is 
not equivalent to a classical spin. The ground state for small 
Ji is a degenerate spin S multiplet, which is not equivalent to 
a single spin-polarized state as predicted by a classical spin 
model (the symmetries of the state are different). Another 
important observation is that, quite generally, in quantum im- 
purity models one needs A ~ Tk to observe sub-gap excita- 
tions well inside the gap (i.e., not at the very edge), as shown 
in Table U Thus in the situations where the spin-flip scatter- 
ing processes play no role, we also do not expect to observe 
any sub-gap peaks deep in the gap. This is in contradiction 
to the results of the classical spin models, which on the one 
hand presume the irrelevance of the Kondo effect due to small 
exchange coupling, yet also predict excitations well inside the 
gap. Such a discrepancy exists even for the relatively large 
spin ^imp = 5/2, which still cannot be considered as a classi- 
cal spin. 

IV. TRANSITIONS IN THE ANISOTROPIC CASE 

In the anisotropic case, the ground-state multiplet with 
S > 1 splits: for axial D < anisotropy the new ground 
state consists of states with the maximal S z — ±5, while for 
planar D > anisotropy the ground state is S z = for in- 
teger S and S z = ±1/2 for half-integer S. The transition 
point at A c in the isotropic models is extended into transi- 
tion lines in the (D, A) plane. For a given value of D, there 
is some gap A where the system makes a transition from a 
"low- 1 S z |" to a "high- 1 S z |" regime, which are the equivalents 
of the Kondo- screened and free-spin regime, respectively. The 
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Figure 2: (Color online) Many-body ground state as a function of 
the gap A and the magnetic anisotropy D at fixed exchange coupling 
J for a range of the impurity spin Si mp in the single-channel case. 
Magnetic field is zero, B = 0. The value A c (circle) corresponds 
to the transition point between the Kondo screened and unscreened 
impurity moment in the isotropic case. 



results of extensive calculations for the single-channel prob- 
lems are summarized in Fig. [2 in the form of phase diagrams, 
while the actual numerical results (including also all sub-gap 
excited states) are shown in the Appendix. 
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Figure 3: (Color online) Left panels: Ground state and sub-gap 
many-body excited states as a function of the magnetic anisotropy 
D for different values of the superconducting gap A for the single- 
channel spin-1 Kondo models. Zero magnetic field, B = 0. Top 
to bottom: "free-spin" regime, transition regime, Kondo screened 
regime. The energies E are plotted relative to the ground-state en- 
ergy Egs and they are rescaled in units of the gap, thus the ground 
state always lies atop the horizontal axis at E — Egs and the contin- 
uum of excitations starts at (E — Eqs) / A = 1. Right panels: En- 
ergies and weights of the sub-gap spectral peaks in the impurity spec- 
tral function. Darker shade corresponds to higher spectral weight of 
the delta peak. 

Let us consider the results for S{ mp = 1 more closely. In 
the left panels in Fig. [3] we plot the ground state as well as the 
sub-gap excited- state energies as a function of the anisotropy 
parameter D for three values of A that correspond in the 



isotropic limit to the Kondo regime, the transition regime, and 
the "free-spin" regime, respectively. The spin-1 ground state 
and excited state (ES) spin multiplets split in the presence of 
the magnetic anisotropy. The degree of splitting is not the 
same as for free multiplets, but rather depends on A and Tk- 
For example, at A/Tk = 8.6 we find 



d(E Sz =i-E Sz =o)/dD = (l\S ] 
« 0.867 



2 

imp, z 



1)- 

0.266 i 



^i 2 mp,,|0) 

0.6, 



(2) 



rather than the free-impurity result 1. The anisotropy effects 
are thus significantly renormalized by the exchange coupling 
of the impurity with the host. This observation is important 
for the interpretation of possible experimental results: from 
the sub-gap excitation spectra, one cannot directly obtain the 
"bare" anisotropy parameters that appear in the Hamiltonian, 
only their effective "renormalized" values. 

The transitions between the ground state and the excited 
states correspond to discrete (delta) sub-gap peaks in the dif- 
ferential conductance spectra. We plot the spectra of these 
peaks in the right panels of Fig. [3j Only the transitions be- 
tween the ground state and the excited states with AS Z = 
±1/2 are observable spectroscopically. Multiple sub-gap 
peaks may be observed, for example, in the "Kondo regime" 
with S z = ±1/2 ground state for D ^ 0. A characteristic 
feature is that some sub-gap peaks may disappear abruptly as 
a function of D/A when the ground state changes. We also 
note that in an interacting superconducting system, it is crucial 
to distinguish between the many-particle states (ground state 
and excited states) and the peaks in the spectral functions as- 
sociated with transition between said states. Namely, to each 
many-particle excitation with energy E, such that AS Z = 
±1/2 with respect to the ground state, corresponds not one 
but two spectral peaks in the single-particle spectral function. 
They are located symmetrically at a; = ±(E — Eqs)', it is 
possible, for example, to end up in exactly the same many- 
particle state by either adding an electron (uj > peak) or by 
removing it (uj < peak). The lack of the distinction between 
the many-particle states and the spectral peaks has led to some 
confusion in the literature. The difference is particularly im- 
portant for interacting systems, where the many-particle states 
cannot always be decomposed into products of single-particle 
levels (quasiparticles). A notable example is the two-channel 
Simp = 1/2 model. 

In fully general multi-channel problems with non-equal J^, 

there are multiple stages of the Kondo screening with differ- 
ed 

ent Kondo temperatures, , i = 1, . . ., TV. Depending on 

the relation between A and all , the system may end up in 
different ground states. As an illustration, in Fig. 01 we depict 
the possible ground states for isotropic N = 2 problems by 
fixing A and plotting the phase diagrams in the ( J\ , J 2) plane. 
The case of = 1/2 is special due to the overscreening 
effects. The second case with S\ 



imp 



1 shows the generic 
phase diagram for all ^imp > 1 and N = 2: in the vicin- 
ity of the equal-coupling line, the ground-state spin changes 
by A^/2 = 1, while for general Ji the multiple Kondo scales 

(i) 

T^ J result in intermediate regimes with only partial impurity 
screening. This is a new feature that is particular to multi- 
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Figure 4: (Color online) Schematic diagram representing the many- 
body ground state as a function of the exchange coupling constants 
Ji and J 2 for the isotropic (D = 0) two-channel S mp = 1/2 and 
Simp = 1 Kondo models at fixed superconducting gap A = 10 -5 . 
Zero magnetic field, B = 0. 

channel problems. In the presence of anisotropy D, the phase 
diagrams in Fig. [H may be extended into the third dimension; 
the GS multiplets with S > 1 split according to the sign of 
D, and effects similar to those represented in Fig. [2 are ob- 
served. (See Fig. [161 for some results of the calculations at 
finite magnetic anisotropy D.) 



V SUB-GAP EXCITED STATES AND SPECTRAL PEAKS 
IN THE MULTI-CHANNEL CASE 



In the presence of additional scattering channels, further sub- 
gap excitations may appear. For N = 2, we find, however, 
that if J 2 is much lower than J\ , the additional excited state 
merges with the continuum and is not observable. For mod- 
erate J 2 1 J\ ratio, the channel-2 excitation peak is still inside 
the gap, but its spectral weight tends to be much smaller than 
that in the dominant i = 1 channel. 

On more symmetric surfaces (adsorption sites), some of 
the exchange coupling constants Ji can be equal, for exam- 
ple pair-wise equal in the presence of a four-fold symmetry 
axis. In this case, the excitations associated with those Ji that 
are equal will be degenerate, thus multiple sub-gap peaks are 
again not expected. 

These results suggest that it is not very likely to observe 
multiple peaks due to coupling to different scattering chan- 
nels. Given that multiple peaks are nonetheless commonly 
observed in experiments 19-21 , we propose that a very likely in- 
terpretation involves the presence of the magnetic anisotropy 
splitting of the sub-gap excitations. An example for TV = 2 
is shown in Fig. [5j For D = 0, the ground state is ^ = 1/2 
and there is an S = 1 excited state associated with channel 
1 and an S = excited state associated with channel 2. For 
D 7^ 0, the triplet excited state splits. In the presence of longi- 
tudinal anisotropy, the peaks will split at most in two, but ad- 
ditional splitting may be induced by the transverse anisotropy 
E(S% — Sy), which is also known to be present in adsorbed 
magnetic impurities 25 . 
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Figure 5: Energies and weights of the sub-gap peaks for the two- 
channel Si m p = 1 Kondo model as a function of the magnetic 
anisotropy. Parameters are A = 10" 5 , Ji = 0.2, and J 2 = 0.1. 
Darker shade corresponds to higher spectral weight. The labels i = 1 
and i = 2 indicate spectral peaks associated with the i-th scattering 
channel. Zero magnetic field, B = 0. 

For real transition-metal impurities on a surface, the ex- 
change coupling constants depend on the energies and the hy- 
bridization constants of the d orbitals. Assuming an adsorp- 
tion site of low symmetry, it is generally more likely that one 
of the exchange constants, say Ji, will be dominant. In strict 
single-channel problems, there is always at least one excited 
state inside the gap. It tends to be located at the gap edge for 
Tk <C A and Tk ^> A, but it is found well inside the gap 
when Tk and A are roughly of the same order of magnitude. 



VI. BEHAVIOR IN THE EXTERNAL MAGNETIC FIELD 

Due to the strongly enhanced spin-orbit interaction on sur- 
faces, the interpretation of multiple peaks in terms of mag- 
netic anisotropy splitting appears very plausible. To experi- 
mentally distinguish between the different possible origins of 
the multiple-peak sub-gap structures in a conclusive way, we 
propose to study the Zeeman splitting of the sub-gap peaks by 
weak magnetic fields (weak enough so that the superconduc- 
tivity is not significantly suppressed; see also Refs. [82l - [85b . 
This is possible because magnetic atoms adsorbed on the sur- 
face of a superconductor are not fully shielded by the Meiss- 
ner effect. Ultra-low-temperature STM's equipped with dilu- 
tion refrigerators are likely to achieve sufficient energy reso- 
lution. The magnetic field splits pairs of the sub-gap states 
with the same | S z | . Since the only observable transitions are 
those with AS Z = ±1/2, each spectral peak may split at 
most in two. In the presence of transverse anisotropy and/or 
transverse magnetic field, the spin symmetry is fully lifted and 
even more complex spectra of excitation peaks can arise. We 
note, however, that if the ground state is spin degenerate (i.e., 
for partially screened impurity), even a small magnetic field 
can fully polarize the residual impurity spin, thus some of 
the sub-gap peaks might not be spectroscopically observable 
since their weights are essentially zero (at T = 0). Some ob- 
served peaks will thus merely shift, rather than split. This be- 
havior is represented schematically in Fig. [6] for the example 
of an ^imp = 1/2 impurity in the single-channel case, while 
the results of a corresponding NRG calculation are shown in 
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Figure 6: (Color online) Schematic representation of the relation be- 
tween the many-body energy levels (shown on the left; Egs is the 
ground-state energy, Ees are the excited-state energies, A is the gap 
beyond which starts the continuum of the conduction-band states) 
and the spectral peaks which can be measured by the STM (shown 
on the right; V > indicates tunneling into "empty" states, i.e., 
the electron addition transitions, while V < corresponds to the 
removal of electrons from the system). Only transitions where S z 
changes by 1 /2 are spectroscopically visible at very low tempera- 
tures. The arrows indicate the spin of the electron being added (pos- 
itive V) or removed (negative V) from the system at the given reso- 
nance. In the absence of the field, Tk > A and Tk < A cases can- 
not be easily distinguished spectroscopically since they both exhibit 
a pair of excitation peaks at uj = ±(Ees — Egs)- In the magnetic 
field, we expect to observe peak splitting for Tk > A, while for 
Tk < A we expect merely a shift of the peak pair concomitant with 
the weight reduction if the experimental temperature is much lower 
than the ^es,i — ^gs scale. The example depicted in the plot cor- 
responds to an Simp = 1/2 magnetic impurity in the single-channel 
case. 



The single-channel S- irnp = 1 model is studied in 
Fig- Eb>c,d) for different values of the magnetic anisotropy 
D. In the absence of the anisotropy, D = 0, the behavior is 
similar to that in the spin-1/2 case: for A/Tk > ("free- 
spin" regime), the peaks merely shift, while for A/Tk < 
(Kondo regime), there is a splitting of the peaks, since the 
transitions are possible from the low-lying doublet state (here 
S z = -1/2) to both S z = -1 and S z = excited states. It 
is worth emphasizing that in the latter case, the peak weights 
are different for the two transitions. A simplified, but intu- 
itive picture is the following: the ground state is a Kondo state 
I Simp, z = — 1 5 T) 5 that is, a bound state of an S z = —1 im- 
purity state and a spin-up conduction-band electron with total 



S z = —1/2. The first excited state has S z = — 1 and can be 
obtained in the process of adding a spin-down electron. The 
second excited state has S z = and it can be reached by 
adding a spin-up electron. Since there is already a spin-up 
electron present in the Kondo ground state, adding a further 
spin-up electron will have a reduced weight as compared to 
adding a spin-down electron. This behavior persists in the 
presence of magnetic anisotropy: the peak corresponding to a 
transition to the S z = state has lower weight than the peak 
corresponding to the transition to the S z = — 1 state in all 
three cases of D = 0, D > 0, and D < 0; see Fig.[7J),c,d). 



S: =1,D=0 




B/A 



Figure 7: (Color online) Energies and weights of the sub-gap spec- 
tral peaks as a function of the external magnetic field for the single- 
channel spin-1/2 and spin-1 Kondo models with different magnetic 
anisotropy terms. 

In the "free-spin" regime, the evolution of the spectral 
peaks in the magnetic field depends on the sign of the mag- 
netic anisotropy; compare the upper panels in Fig. [7}),c,d). 
The isotropic D = and easy-axis D < cases are similar: in 
the presence of the field, the S z = — 1 state will be the ground 
state and the transitions are only possible to the S z = —1/2 
excited state, thus a single sub-gap peak pair is observed. The 
easy-plane D > case is more interesting. Now the ground 
state (for small magnetic fields) is S z = 0, thus the transitions 
to both S z = +1/2 and S z = —1/2 excited states are possi- 
ble; we thus see two pairs of sub-gap spectral peaks, that is, a 
total of four sharp peaks. For very large Zeeman splitting, the 
S z = —1/2 excited state will become the new ground state 
of the system; at this point only the transition to the 5^=0 
sub-gap will be possible and a single pair of peaks will remain 
in the impurity spectral function. 

In Fig.[8j we finally plot the field dependence of the spectral 
peaks of a two-channel spin-1 impurity for different magnetic 
anisotropies D. In all three cases, the ground state for B > 
is non-degenerate, S z = —1/2, and one may only observe 
the transitions to S z = —1 and S z = excited states. As 
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evidenced by the results in Fig. [8l the sign and magnitude of 
D can be easily determined from the shifts. The peak due 
to the second weakly coupled conduction channel is always 
weaker as compared to the peaks associated with the dominant 
screening channel. 
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Figure 8: Energies and weights of the sub-gap spectral peaks as a 
function of the magnetic field B for the two-channel S mp = 1 im- 
purity. Parameters are as in Fig. \5\ The spectral peaks are fully 
spin-polarized, as indicated by the arrows. 



VII. EFFECTS OF THE INTER-IMPURITY EXCHANGE 
INTERACTION 

When adsorbed magnetic impurities are brought to- 
gether, for example by controlled manipulation using the 
tip of the STM, their mutual interaction will change 
the LDOS signatures measured by the scanning tunneling 
spectroscopy^*^ - — . The theory of the inter-impurity interac- 
tions in the normal case has been worked out using both sim- 
plified model Hamiltonians and ab initio calculations 9 ^ - —. In 
the superconducting case, the calculations of the sub-gap exci- 
tation spectra in impurity dimers have been performed mainly 
for classical impurity spin a 96 ! 97 . Recently, some calculations 
for coupled quantum impurities have been performed in the 
context of double quantum dots 9 ^. Here we study the problem 
of two quantum impurities, each coupled to a separate con- 
duction band with a superconducting gap (with equal A), and 
interacting via an isotropic Heisenberg Hamiltonian: 



JS 



imp,l " ^imp,2- 



(3) 



We first discuss the simplest case of two 5i mp = 1/2 impu- 
rities. The numerical results for the sub-gap states are shown 
in Fig. [9] The parameters are chosen such that in the absence 
of coupling, each impurity is the Kondo screened regime with 
an 5 = ground state and one 5 = 1/2 sub-gap excited 
state (two-fold degenerate due to spin); for the two-impurity 
system, the ground state is thus a singlet (we will call it the 
"Kondo singlet"), and there are two degenerate 5 = 1/2 ex- 
cited states (four-fold total degeneracy). As the antiferromag- 
netic ( J > 0) exchange coupling is turned on, a new singlet 
state emerges from the continuum; this state can be interpreted 







~S=l/2 




■ "S=l 





< 0.6 





00060 

4xl0" 5 8x10"^ 



A/T K =0.86 



S. =1/2 

imp 

-4- 



0.0005 



0.001 



Figure 9: (Color online) Ground state and sub-gap many-body ex- 
cited states in the two-impurity problem with 5i mp = 1/2 impurities 
coupled by Heisenberg exchange coupling J. The inset is a close-up 
view of the region where the two different singlet ground states cross: 
the arrow shows the evolution of the local AFM-ordered spin- singlet 
ground state as J is increased. The gap is A = 10 -5 . 



as arising from the antiferromagnetically ordered local singlet 
state formed by the two impurity spins, and we will there- 
fore denote it as the "AFM singlet". It should be emphasized 
that this state is not directly related to the 5 = 1/2 excited 
states (which are actually present in the sub-gap spectrum at 
the same time as the new AFM singlet); the AFM singlet state 
should thus not be interpreted as emerging from the coupling 
of the 5 = 1/2 sub-gap localized states, but rather as aris- 
ing from the local inter-impurity singlet state. We further- 
more emphasize that the two 5 = 1/2 excited states remain 
degenerate (this holds, in fact, for all values of J). As J is 
increased beyond J « 6 x 10 -5 , the Kondo and AFM sin- 
glet states cross and exchange their roles as the ground and 
the excited state. At still higher J, both the excited Kondo 
singlet state and the two 5 = 1/2 excited states merge with 
the continuum and are no longer observable. For large anti- 
ferromagnetic inter-impurity coupling, the dimer behaves as a 
non-magnetic object and therefore does not have any excita- 
tions deep inside the gap. Such behavior seems to be present 
for Cr dimers (in the configuration "Cr dimer II" with small 
inter- atom separation) 19 . 

Equally interesting is the case of ferromagnetic Heisenberg 
coupling, also shown in Fig.O In this case, a new "FM triplet" 
sub-gap state emerges in the sub-gap spectrum. This state de- 
creases in energy until it replaces the Kondo singlet as the new 
ground state. Two important observation can be made: (i) the 
level crossing between the FM triplet and the Kondo singlet 
occurs for a much larger (by an order of magnitude) absolute 
value of the Heisenberg coupling as the level crossing between 
the AFM singlet and the Kondo singlet; (ii) with increasing 
I J\, the Kondo singlet and the degenerate 5 = 1/2 excited 
states evolve only slowly and remain deep inside the gap even 
for very large ferromagnetic exchange coupling. This sug- 
gests that in the case of ferromagnetic dimers, we are more 
likely to observe some sub-gap spectral peaks. It is possible 
that the Mn dimers (even in the configuration "Mn dimer II" 
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with small inter-atom separation^) have ferromagnetic spin 
coupling. 




D/A=-l 



0+ 

-0.0006 



-0.0004 



# ■ — 1 — ■ — * m )ood)ooocb — 
-0.0002 0.0002 



Figure 10: (Color online) Ground state and sub-gap many-body ex- 
cited states in the two-impurity problem with 5i mp — 1 impurities 
coupled by Heisenberg exchange coupling J. The gap is A = 10 -5 . 

We now turn to the case of two isotropic 5i mp = 1 impu- 
rities, each coupled to a single channel; see the upper panel 
in Fig. |T0j We again consider the parameter regime where 
each impurity is Kondo-screened. This time, however, the im- 
purities have residual spin 1/2. In the absence of the inter- 
impurity coupling, the ground state is degenerate and com- 
posed of a singlet and a triplet "Kondo state", and there are 
degenerate excited states with 5=1/2 and 5 = 3/2. The 
ground- state degeneracy is lifted by a small J, thus for small 
J the ground state is either a singlet (for J > 0) or a triplet 
(for J < 0); this splitting occurs on the energy scale of almost 
"bare" J; see the close-up view in Fig.Qj] For antiferromag- 
netic exchange coupling, there is an additional singlet state 
that arises from the local inter-impurity singlet state in which 
two 5 = 1 spins are rigidly antiferromagnetically ordered. 
This state is different from the singlet Kondo state, which 
arises from AFM ordering between two residual 5 = 1/2 
extended objects. This is another confirmation of the nature 
of the sub-gap singlet states, which we had already discussed 
in the case of 5i mp = 1/2. For sufficiently large AFM cou- 
pling, the singlet states cross, and for very large coupling we 
again find that there are no sub-gap excited states. Note also 
the similarity between the J > behavior in 5i mp = 1/2 and 



5im P = 1 models. 

For ferromagnetic coupling, the situation is again similar to 
what we had observed for 5i mp = 1/2 impurities: the fer- 
romagnetically ordered 5 = 2 local object becomes the new 
ground state only for very large Heisenberg coupling. The 
transition occurs at J w —6 x 10 -5 , that is, at essentially 
the same value as for 5i mp = 1/2. In addition, we again ob- 
serve that some sub-gap excited states are typically found in 
the sub-gap spectrum even for large \J\. 
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Figure 11: (Color online) Close-up view of the small- J region in 
Fig.[[0l 

Finally, we discuss how the sub-gap excitation spectra 
change in the presence of the magnetic anisotropy. In the two 
lower panels in Fig. [lOl we show the results for easy-plane 
and easy-axis magnetic anisotropy for the 5i mp = 1 impuri- 
ties. The results are in general similar to those for the isotropic 
model, but there are some notable differences. 

For AFM coupling , we find that the ground state is an S z = 
state (i.e. , the "singlet" state) and that there are no sub-gap 
excited states for large enough J. In fact, it may be noticed 
that the results for J > are remarkably similar irrespective 
of the sign and strength of the magnetic anisotropy. 

For small ferromagnetic exchange coupling, the 5 = 1 
ground state splits: for easy-plane anisotropy the ground state 
is S z = 0, while for easy-axis anisotropy it is S z = ±1. 
Thus the combined inter-impurity states can be understood 
within a two- stage splitting model: in the first step, we con- 
sider how the states combine due to the Heisenberg exchange 
coupling; in the second step, we consider how the resulting 
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states split due to the magnetic anisotropy. Such separation is 
possible because in the parameter regime under consideration, 
the scale of the exchange splitting (up to several times 10 ~ 4 ) 
is larger than the scale of the magnetic anisotropy (fixed at 
\D\ = 10 -5 ). For strong FM coupling, the result depends on 
the type of anisotropy: the ground state arises from an S = 2 
state, but for the easy-plane case the actual ground state is a 
non-degenerate S z = state, while for the easy-axis case the 
ground state is a twofold-degenerate S z = ±2 state. The tran- 
sition between the small- \J\ and large- \J\ ground state occurs 
at a value of J that depends on the magnetic anisotropy. 

VIII. CONCLUSION 

We have studied single magnetic impurities on supercon- 
ductor surfaces as well as their dimers. We have emphasized 
the importance of using quantum impurity models to describe 
the sub-gap excitation spectrum of the magnetic adatoms. The 
isotropic models exhibit a transition between ground states 
with different degrees of Kondo screening of the impurity 
spin; this depends on the values of the Kondo exchange cou- 
pling constants and the ratios between the resulting Kondo 
temperatures and the BCS superconducting gap. The sub-gap 
excited states only appear in cases where at least one of the 
Kondo temperature scales is of the same order as the super- 
conducting gap. If all coupling constants are small (i.e., if the 
Kondo effect is absent), there are no sub-gap excited states 
deep inside the gap, which is at odds with the predictions 
from the classical spin models. The ground state and sub- 
gap excited states with S > 1 are split in the presence of the 
magnetic anisotropy. It is found that the splitting is strongly 
renormalized by the Kondo screening, thus the model needs 
to be studied using non-perturbative techniques, such as the 
numerical renormalization group. We find that the weight of 
the sub-gap peaks is the largest for strongly coupled Kondo 
channels, but becomes lower for weakly coupled channels. 
In addition, the sub-gap excitations associated with weakly 
coupled channels tend to appear near the gap edge, and may 



therefore be difficult to observe. This finding suggests that 
multiple-peak sub-gap excitations likely arise from the inter- 
nal structure of the impurity (spin-orbit coupling leading to 
the magnetic anisotropy). We have also explored the inter- 
dot exchange coupling for impurity dimers: for strong anti- 
ferromagnetic coupling, no sub-gap excitations are present, 
while even for relatively strong ferromagnetic coupling some 
sub-gap peaks may be observable. We have shown that the 
external magnetic field has a sizable effect on the sub-gap ex- 
citation spectra of impurities. Exploring the magnetic prop- 
erties of impurities using the field dependence of the sub-gap 
peaks constitutes a worthwhile experimental challenge. Fur- 
thermore, using a break-junction setup with superconducting 
contacts should make it possible to map the excitation spectra 
in the entire (D,B) plane at the same time, providing further 
means to test the predictions of this work. 
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Appendix A: Sub-gap excitation spectra 

This appendix contains additional figures detailing the re- 
sults discussed in the main text. Figures [T2l and [T3l show the 
sub-gap states for one-channel and two-channel Kondo mod- 
els from which the transition points between the (partially) 
screened and unscreened regimes can be extracted; see Ta- 
ble [Band Figure [2 in the main text. Fig. ITU shows the effect 
of the magnetic anisotropy on the sub-gap states, that is, the 
splitting of the degenerate spin multiplets. Figures [T5l and [T6l 
contain the results for the two-channel S = 1/2 and 5 = 1 
Kondo model in the (Ji, J2) plane, which serve to establish 
the phase diagram in Fig. 01 
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Figure 13: (Color online) Ground- state and sub-gap many-body excited states as a function of the superconducting gap A for two-channel 
(N — 2) Kondo models with different impurity spin S mp . The NRG calculations have been performed with the discretization parameter 
A = 4. The results given in Table 1 for the two-channel case have been calculated by performing further calculations using a twist parameter 
z — 0.5 and averaging the results for z — 1 and z — 0.5. The tabulated results are then accurate within a few percent for both the single- 
channel and the two-channel case. Note the presence of the additional S — S'imp — 1/2 excited-state multiplet inside the gap. For non-equal 
Ji, this excited state can become the ground state in some parameter regimes; see Fig. 01 




Figure 14: (Color online) Ground- state and sub-gap excited states as a function of the magnetic anisotropy D for different values of the 
superconducting gap A for the single-channel spin-1, spin-3/2, and spin-2 Kondo models. The results of these (and similar) calculations have 
been used to establish the schematic phase diagrams shown in Fig. [2 
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Figure 15: (Color online) Ground state and sub-gap excitations as a function of the parameter Jo for the two-channel Simp = 1/2 Kondo 
model. We define Ji = Jo cos and J 2 — Jo sin 0; the angles of the direction in the ( Ji , J 2) plane are given (in degrees) as the titles of the 
subfigures. The gap is A = 10 -5 . The results of these (and similar) calculations have been used to establish the schematic phase diagram for 
the S'imp = 1/2 case in Fig. [4] of the main text. The two-channel calculations with non-equal Ji have to be performed with a higher value of 
the NRG discretization parameter A = 8, otherwise they are intractable. No twist-parameter averaging has been performed here; nevertheless, 
the results are still qualitatively correct. 
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